The Instanton/Sphaleron Mechanism of Prompt Gluon Production 
in High Energy Heavy Ion Collisions at RHIC 
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We argue that if the growing part of hadron-hadron cross section (described phenomenologically by 
the supercritical soft Pomeron) is due to instanton/sphaleron mechanism, one should find certain 
qualitative features of the produced cluster which differ from the usual string fragmentation. Fur- 
thermore, we suggest that this mechanism should be even more important for heavy ion collisions 
in the RHIC energy domain. Large number of parton-parton collisions should result in hundreds 
of produced sphaleron-like gluomagnetic clusters per unit rapidity. Unlike perturbative gluons (or 
mini-jets), these classically unstable objects promptly decay into several gluons and quarks in mini- 
explosions, leading to very rapid entropy generation. This may help to explain why the QGP seem 
to be produced at RHIC so early. We further argue that this mechanism cannot be important at 
higher energies (LHC), where perturbative description should apply. 



l.At high energies s > 10 3 GeU 2 hadronic cross sec- 
tions (pp, pp, np, Kp, jN and even 77) slowly grow with 
the collision energy s. This behavior can be parame- 
terised by a Regge pole, the so called soft Pomeron (see 
e.g. [GJ). In this Letter we will not address very high s 
and therefore use only the logarithmic fit 



o-hh'(s) = o- h h'(s ) + log(s / so)X hh , A + 



(1) 



ignoring both the question of whether it is indeed a Regge 
pole, as well as other Reggions leading to contribution 
decreasing with energy. For estimates below we use pa- 
rameters from the latest Particle Data Group fits 0, 
which give the "Pomeron intercept" and a constant equal 
to pp,pp collisions, A = a(0) — 1 = 0.093(2), Xnn — 
18.951(27) mb. 

A qualitative difference between constant and logarith- 
mically growing parts of the cross section will be em- 
phasised. The former can be explained by prompt color 
exchanges, as suggested by Low and Nussinov Q long 
ago. The growing part of the cross section cannot be 
generated by t-channel vector exchanges and is associated 
with prompt (prior to formation of strings) production of 
some objects, with log(s) coming from longitudinal phase 
space. Perturbative QCD describes gluon production, by 
processes like the one shown in Fig.[j](a), which can be it- 
erated in the t-channel in ladder-type fashion resulting in 
a BFKL pole M. Although its intercept is much larger 
than A mentioned, it is consistent with much stronger 
growth seen in hard processes at HERA: thus it is often 
called a "hard pomeron" . 

The physical origin of growing cross section remains 
an outstanding open problem: neither the perturbative 
resummations nor existing non-perturbative models are 
really quantitative. It is hardly surprising, since the scale 
at which soft Pomeron operates (as seen e.g. from the 
Pomeron slope a'(0) w 1/(2 GeV) 2 ) is the semi-hard or 
"substructure scale" Q 2 ~ 1 — 2GeU 2 , which is notori- 
ously difficult for theorists because it is simultaneously 
the lower boundary of pQCD (serving therefore as the 



cut-off Pmin already mentioned), as well as the upper 
boundary of low energy effective approaches like chiral 
Lagrangians or Nambu-Jona-Lasinio model. At the same 
time, a number of objects/phenomena are naturally as- 
cribed to this scale: "constituent quarks" , flux tubes 
(or QCD strings) and their junctions, to name a few. 
We do not have a quantitative description of flux tubes 
(other than lattice QCD), but constituent quarks and 
related issues can be well understood in an instanton liq- 
uid model, see review || . Its primary parameters are the 
number density of instantons (plus anti-instantons) and 
their average size, determined long ago [Q, from QCD 
phenomenology to be n « 1 fm and small average size 
of p k, 1/3 fm leading to vacuum diluteness np 4 ~ 10~ 2 . 
Amazingly, with those two numbers one can get truly 
quantitative description of correlators, form-factors and 
other hadronic parameters^. 

Application of the instanton- induced dynamics to high 
energy hadronic collisions have been suggested recently 
P,^0|. One important precursor has been the Kharzeev- 
Levin work [|| in which contribution to A of scalar color- 
less states - the sigma meson and the scalar glueball - has 
been non-perturbatively evaluated. Two last works in 
fiof have benefited from deep insights obtained a decade 
ago in studies of instanton-induced processes in elec- 
troweak theory, see Jll| and references therein. In these 
works the growing part of the hh cross sections is due to 
prompt multi-gluon production via instantons, or more 
accurately, via colored gluonic clusters called sphalerons, 
see Fig.|(b). 

Among qualitative features of this theory is the ex- 
planation of why no odderon appears (instantons are 
SU(2) objects, in which quarks and antiquarks are not 



*For recent example see || where vector and axial correlators 
obtained from the r decays are explained literally within their 
error bars, or withing few percent accuracy. 
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really distinct), an explanation of the small power A (it 
is proportional to "instanton diluteness parameter" np A 
mentioned above), the small size of the soft Pomeron 
(governed simply by small size of instantons mentioned 
above, p ~ 1/3 fm). Although instanton- induced ampli- 
tudes are proportional to small "diluteness" factor, there 
is no extra penalty for production of new gluons: thus 
one should expect instanton effects to beat perturbative 
amplitudes of sufficiently high order. This generic idea is 
also behind the present work, dealing with prompt multi- 
gluon production. 




FIG. 1. (a) A typical inelastic perturbative process: 
two t-channel gluons collide, producing a pair of gluons; (b) 
Instanton-induced inelastic process incorporate collisions of 
multiple t-channel gluons with the instanton (the shaded cir- 
cle), resulting in multi-gluon production. The intermediate 
stage of the process, indicated by the horizontal dashed lines, 
corresponds to a time when outgoing glue is in the form of 
coherent field configuration - the sphaleron. 

Technical description of the process is split into two 
stages. The first (at which one evaluates the cross sec- 
tion) is the motion under the barrier, described by Eu- 
clidean paths approximated by instantons. Their in- 
teraction with the high energy colliding partons results 
in some energy deposition and subsequent real motion 
above the barrier. At this second stage the action is 
real, and the factor \exp(iS)\ — 1, so it does not af- 
fect the cross section and is need only to detail the final 
state. The relevant Minkowski paths start with configu- 
rations close to QCD analogs of electroweak sphalerons 
p2| , static spherically-symmetric clusters of gluomag- 
netic fieloj^]. Their mass in QCD is also determined by 
the isnatnton size 



M sph w - ~ 2.5 Gel/ (2) 

9\P)P 

Since those field configurations are close to classically 
unstable saddle point at the top of the barrier, they roll 
downhill and develop gluoelectric fields. When both be- 
come weak enough, solution can be decomposed into per- 
turbative gluons. This part of the process can also be 
studied directly from classical Yang-Mills equation: for 
electroweak sphalerons it has been done in Refs Jl3| , cal- 
culation for its QCD version is in progress |Q. While 
rolling, the configurations tend to forget the initial im- 
perfections (such as a non-spherical shapes) since there is 
only one basic instability path downward: so the result- 
ing fields should be nearly perfect spherical expanding 
shells. Electroweak sphalerons decay into approximately 
51 W,Z,H quanta, of which only about 10% are Higgses, 
which carry only 4% of energy. Ignoring those, one can 
make crude tentative estimate of mean gluon multiplicity 
per sphaleron decay, by simple re-scaling of the coupling 
constants 

a 2 

< N q >«< N W z> yelectr ™ e " k ~ 3 - 4 (3) 
9qcd(p) 2 

The spectrum (also derived from a solution [^3|,0) has 
a wide maximum and can roughly be approximated by 
thermal one, with a temperature of about Tq ~ 300M eV, 
see @. 

2. The first points we would like to make in this Letter 
are some suggestions of how one can experimentally test 
this scenario, by some qualitative effects. 

Note that the fate of the produced sphalerons is differ- 
ent in hh and A A collisions: in the former case they decay 
in the confining vacuum, in the latter into a deconfined 
media (see below). Some fraction of produced clusters 
have net zero color and can directly form glueballs, with 
J p = + ,0~ ^. The scalar isoscalar channel has been 
considered first in ||: however, it can only account for 
a fraction^] of prompt production. Most of the promptly 
produced gluon clusters have non-zero net color, and the 
thus have to be connected by color flux tubes to other 
partons. 

This clearly makes their direct observation difficult, 
but not hopeless: we briefly describe two particular ideas 
of how it can be done. The scalar glueball candidate 
/o(1700) decay into rjrj,KK and only a little bit into 
7T7T. We do not yet have experimental pseudo-scalar glue- 
ball candidate, while lattice predicts it to be right at 
the mass of the sphaleron (2). However, as noticed by 



2 Those can be obtained from known electroweak solutions 
in the limit of infinitely large Higgs self-coupling. 



3 But not a channel with e.g. 2 + quantum numbers, which 
does not classically couple to instantons . 

4 In terms of pomeron intersept, it is A + w 0.05 for scalar 
glueball and sigma together || , while (including shadowing) 
the experimental total is « 0.16 ]33[. 
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Bjorken |^T|, the n c decay has three distinct 3- meson 
modes, KKn, nmr, ti'tttt, with about 0.05 branching each: 
those fit well to the idea that they come directly from 
the 't Hooft instanton- induced Lagrangian uuddss. Pre- 
sumably the instanton-induced decays modes of the 0~ 
glucball should prefer the same 3 channels. 

One may speculate further, and suggest that scalar 
(pseudo-scalar) projections of the sphaleron may still fol- 
low the same scalar (pseudo-scalar) glueball decay pat- 
tern, even while the total color is non-zero. The pat- 
tern of enhanced production of if, 77, K via strange part 
of 't Hooft Lagrangian leads to a specific feature of the 
final state. Indeed, when r)',r],Ks decay into 5,3,2 pi- 
ons, respectively, all of them are produced much later 
than the average pion production time. They are dif- 
ferent from others in one important aspect: they do 
not participate in Bose-Einstein (or HBT) correlations. 
Its strength is traditionally expressed in terms of the so 
called Xhbt = (1 — f) 2 , where / is the fraction of pions 
coming from long-lived^ sources. In minimal bias pp, or 
heavy ion collisions with any multiplicity, or in the e + e~ 
reactions the usual value Xhbt ~ 0.5. However for high 
multiplicity pp collisions experiments show that the in- 
tensity of the correlations decreases substantially, to only 
Xhbt ~ 0.2. As far as we know, this effect has not been 
explaned: see discussion of data and proposed sugges- 
tions in |3^]. Although clearly much more studies are 
needed, it may indicate that promptly produced hadrons 
have an origin other than the usual string fragmentation. 

Another possible approach is based on the (so far ig- 
nored) topological properties of instantons/sphalerons. 
Roughly speaking, each sphaleron has an option to roll 
down in two directions, selecting two possible orienta- 
tions of its gluoelectric field relative to gluomagnetic 
one. Parity is of course conserved in QCD: but on the 
event — by — event basis large fluctuations may appear in 
P- and CP-odd kinematical observables specified in [^8| . 

3. We now turn to heavy ion collisions. Recent ex- 
periments at Relativistic Heavy Ion Collider (RHIC) at 
Brookhaven National Laboratory, taken during its first 
run in summer 2000 and reported recently at Quark Mat- 
ter 2001 conference [Q, have shown that heavy ions col- 
lisions (AA) at highest energies significantly differ both 
from the hh collisions and the AA collisions at lower 
(SPS/AGS) energies. Many features of these data are 
quite consistent with the Quark-Gluon Plasma (QGP) 
(or Little Bang) scenario [Q, in which entropy is pro- 
duced promptly and subsequent expansion is close to adi- 
abatic expansion of equilibrated hot medium. 

Already the very first multiplicity measurements re- 
ported by PHOBOS collaboration |Q have shown that 
particle production per participant nucleon is no longer 



5 Defined relative to h/AE where AE is the energy resolution 
of the detector. 



constant, as was the case at lower (SPS/AGS) energies, 
but grows more rapidly. This behaviour may be due to 
long-anticipated pQCD processes, leading to perturba- 
tive production of new partons. Unlike high p t processes, 
those are (directly undetectable) "mini-jets" . Their pro- 
duction and decay was discussed in Refs JjiJ, and also 
used in widely used event generator HIJING ^0|. Its 
crucial parameter is the cutoff scale p m i n — 1.5 — 2 GeV: 
if fitted from pp data to be, it leads to predicted mini-jet 
multiplicity dN g /dy ~ 200 for central AuAu collisions at 
-y/s = 130 AGeV. If those fragment independently into 
hadrons, and are supplemented by "soft" string-decay 
component, the predicted total multiplicity was found 
to be in good agreement with the first RHIC multiplic- 
ity data. Because partons interact perturbatively, with 
their scattering and radiation being strongly peaked at 
small angles, their equilibration is expected to be rela- 
tively long . 

However, next set of RHIC data reported in |H| have 
provided serious arguments against the mini-jet scenario, 
and point toward quite rapid entropy production rate and 
early QGP formation. 

(i) If most of secondaries come from independent mini- 
jets fragmentation, there would be no collective phenom- 
ena such as transverse flow related with the QGP pres- 
sure. However, such effects are very strong at RHIC. In 
particular, STAR collaboration have observed very ro- 
bust elliptic flow [^2| , which is in perfect agreement with 
predictions of hydrodynamical model p3| assuming equi- 
librated QGP with its full pressure p w e/3 above the 
QCD phase transition. This agreement persists to rather 
peripheral collisions, in which the overlap almond-shaped 
region of two nuclei is only a couple fm thick. STAR and 
PHENIX data on spectra of identified particles, espe- 
cially p,p, indicate spectacular radial expansion, also in 
agreement with hydro calculations p3| . 

(ii) Spectra of hadrons at large pt, especially the 7r° 
spectra from PHENIX, agree well with HIJING for pe- 
ripheral collisions, but show much smaller yields for cen- 
tral ones, with rather different spectra both in shape and 
composition. Moreover, those agree weel with hydro pre- 
dictions which had been established at low p t previously. 
It means that not only long-anticipated "jet quenching" 
is observed, it seems to be as large as it can possibly 
be^]. For that to happen, the outgoing high-p t jets should 
propagate through matter with parton population much 
larger than the abovementioned minijet density predicted 
by pQCD (HIJING). 



Jets originating from the surface outward is very difficult 
to quench, and thus the suppression factor of about ~ 1/10 
is difficult to decrease further, whatever happens in dense 
matter. Counting from expected Cronin effect (which in pA 
collisions is about factor 2 at pt in question), the observed 
suppression is not far from such number. 
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(iii) Curious interplay between collective and jet ef- 
fects have also been studied by STAR collaboration, in 
form of elliptic asymmetry parameter v^ipt). At large 
transverse momenta pt > 2 GeV the data behave accord- 
ing to predictions of jet quenching model pffl , indicating 
gluon multiplicity several times larger than HIJING pre- 
diction. Moreover, the result is in fact consistent with the 
maximal possible value evaluated from the final entropy 
at freeze-out, (dN/dy) w ~ 1000. 

In this Letter we propose a non-perturbative solution 
to this puzzle. But before we come to it, let us also 
mention its alternative: significantly lower cutoff scale 
in excited matter, as compared to p m in = 1-5 — 2 GeV 
fitted from the pp data. It may lead to larger perturba- 
tive cross sections, both due to smaller momenta transfer 
and larger coupling constant. As argued over the years 
(see e.g. one of the talks the QGP is a new phase 
of QCD which is qualitatively different from the QCD 
vacuum: therefore the cut-offs of pQCD may have en- 
tirely different values and be determined by different phe- 
nomena. Furthermore, since QGP is a plasma-like phase 
which screens itself perturbatively |l6[ , one may think 
of a cut-offs to be determined self- consistently from re- 
summation of perturbative effects. These ideas known as 
self- screening or QGP saturation were discussed in Refs. 
||25| . Although the scale in question grows with temper- 
ature or density, just above T c it may actually be smaller 
than the value 1.5-2 GeV observed in vacuum. Its direct 
experimental manifestation may be deformation of dilep- 
ton spectra, wheih can be well described by decreasing 
"duality scale" @. 

4. In order to specify the magnitude of new production 
mechanism one can study dependence on the impact pa- 
rameter b. This dependence of the pseudo-rapidity den- 
sity at mid-rapidity, measured at RHIC [|l5| can be very 
accurately described by simple parameterization p7| 



dn AuAu (n = Q,b)/dr) _ (1 - x) 
driNNiv = 0)/dr) 2 



Npart(b) + xN coU (b) (4) 



where the average number of participants N part (b) and 
NN collisions N co u(b) are calculated in standard Glauber 
model. 

The key is new (b-independent) parameter x(s), defin- 
ing which fraction of NN collisions scales differently in 
AA. We propose to identify x(s) with the growing part 
of the NN cross section discussed above, namely 



x(s) = A XnN log(s/s ) 
o~hh'(so) 



(5) 



Note that two phenomenological values fitted at two 
RHIC energies x{*/s = 56Gev) = 0.05 ± 0.03, x(*Js = 
130 Get;) = 0.09 ± 0.03 @ are both well reproduced if 
one selects the threshold value at sq = 1000 GeV 2 , the 
position of the NN cross section minimum. Further- 
more, because this s is above the highest SPS energy, 
it explains why this component has not been seen be- 
fore. This identification is due to the picture of prompt 



production of some objects - minijets or sphalerons - in 
partonic collisions f\. 

Partons which participate in such interaction should be 
appropriately normalized at the scale y? ~ 1 — 2 GeV 2 . 
Constituent quark models of 60 's would count only them, 
so Nbaryons _ g Using parton densities derived from 
structure functions one finds that (at scale under dis- 
cussion) sea can be neglected, but gluons do not. With 
significant uncertainties, for RHIC energy one can inte- 
grate structure functions for x > 0.01 and get roughly 
additional 3 gluons, leading to sa 6 [| 

The inelastic hh' cross section can be schematically 
written in a simple multiplicative form 



a hh = N^N p l (a° pp + a^logs + 0{log 2 {s)) 



(0) 



For simplicity of presentation, we ignore the difference 
between qq, qq, qg, gg cases, as well as possible depen- 
dence on quark flavor. Here N p is the number of partons 
per hadron and a pp , a pp are the parton-parton cross sec- 
tions, without and with prompt production. In what 
follows we ignore the former and only concentrate on the 
latter part, normalizing it to the observed soft Pomeron 
growth crL = Xnn/(N p ) 2 . Bypassing dynamical cal- 
culations pXj ], we then estimate the probability of the 
sphaleron production directly from data, by assuming 
it to be the dominant process behind the logarithmic 
growth of the cross section. 

It means that in mean parton-parton collision, the 
cross section per rapidity of prompt productionn 



do 



prompt 



dy 



AatatA 



0.005 fm 2 



(7) 



Now we evaluate the total number of parton-parton col- 
lisions in central AA collision. Unlike the total cross 
section, it is not just a multiplicative expression: nuclear 
geometry leads to 



N pp collisionsiAA) ~ (A * A^) 4/3 ~ 10 4 



(8) 



where in numerical estimate wc have used A — 200. As- 



7 The const(s) part of the cross section, which is associated 
with color exchanges, should scale as the number of partic- 
ipants because, no matter how many exchanges took place, 
each outgoing parton pulls out only one color flux tube per 
quark (or 2 per gluon). 

8 Detailed evaluation of semi-hard partonic cross sections 
from (i) the growing part of all hadron-hadron scattering cross 
sections, (ii) elastic amplitudes and (iii) p,7r,7 structure func- 
tions will be reported elsewhere ps) |. 

9 Note a surprisingly small, factor 1/100, compared to geo- 
metric cross section np 2 . In instanton-based theory it origi- 
nates directly from the first power of instanton diluteness 
of such magnitude. 
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suming simple factorization of the cross sections^. Com- 
bining these two simple ingredients we now estimate the 
total density of "promptly-produced objects" (mini-jet 
pairs or sphalerons) in AA collisions per unit rapidity 



dN, 



prompt 



x 



N N l 



(III 



! W 2/3 



200 



(9) 



Presumably one can still treat these objects as produced 
independently, since the number of available cells in the 
transverse plane N ce u s = (R/p) 2 « 400 is still larger than 
this maximal sphaleron number. 

The number of "promptly produced objects" estimated 
above is rather close to mini jet-productionp] calculated 
with HIJING |^(|. Furthermore, multiplying it by the 
transverse energy (2/3)m sp ^, we find that our prompt 
production should result in roughly dE t /dy ~ 400 GeV 
of transverse energy, again comparable to HIJING pre- 
dictions. So a critical reader may ask whether actually 
anything has been gained, by substituting one hypothet- 
ical mechanisms of prompt production - the mini-jet sce- 
nario - by another one, based on instantons/sphalerons. 
Indeed, provided both are similarly normalized to grow- 
ing part of the pp collisions and then scaled to AA case, 
we get about the same number of semi-hard events and 
the sama excitation energy. 

Our first (theoretical) answer is that the suggested sce- 
nario suggests an explanation to the semi-hard scale in- 
volved, derived from the well known vacuum instanton 
parameters, while in pQCD the cutoff should be just 
guessed or fitted. Furthermore, it implies detailed mi- 
croscopic knowledge of the specific gluon field configura- 
tion involved, not just estimate of a number of gluons 
produced. 

The second (pragmatic) answer is that these two sce- 
narios differ significantly in the amount of the entropy 
produced. The minijets are just plane waves: they are 
classically stable and weakly interacting. The sphalerons 



10 Note that we are still very far from unitarity constraints. 
Inside the tube with instanton radius Tvp 2 we find about .67 
partons in a nucleon and 3.6 in Au: so even factorized cross 
section lead to interaction probablity of only 1/200 and 1/10 
respectively, much less than 1. It does not mean however 
that factorization is accurate: we use it only as an estimate. 
For instanton processes presence of extra partons lead to ex- 
tra factors - Wilson lines - in the amplitude, but averaging 
over instanton collective variables (such as color orientations) 
may upset factorization. This question deserves quantitative 
study. Note also, that partons found at the same position in 
transverse plane most likely come from different nucleons, so 
position and color correlations between them are likely to be 
small: so no assumptions about wave functions are probably 
needed here. 

11 As should be expected, since minijet models fit the pp 
cross section as well. 



are unstable, a kind of resonances existing already at clas- 
sical level. They explode into spherical expanding shells 
of strong field, which rapidly sweep the whole volume 
and may convert it into Quark-Gluon Plasma, in which 
the charge is screened rather than confined [M. The 
"initial temperature" of gluons produced from sphaleron 
decay T ~ 300 Mew indicated above is definitely above 
the critical value. Most important, the produced entropy 
is several times larger than for minjets, as recent RHIC 
data seem to indicate. 

In heavy ion collisions at RHIC the QGP is supposed 
to exist at RHIC for several fm/c, much longer than the 
sphaleron lifetime T sp h ~ 1/p. If so, partons produced 
do not hadronize immediately (as for hh collisions) but 
decay into 3-4 gluons, plus 0-6 quarks^ and start real 
equilibration. 

Phenomenologically, comparing dN spha i erons / dy « 
200 sphalerons to dN g i uons /dy « 1000 one sees that 
about 5 partons/sphaleron would produce the right 
amount of entropy, that about 5 partons/sphaleron 
would do the job, which is concievable. In order to test 
the conjectured mechanism experimentally, one may try 
to infer gluon/quark ratio at early time from dilepton 
production. Another possibility is to look at event-by- 
event fluctuations following from clustering at the pro- 
duction stage. 

5. Finally, let us briefly discuss what we predict should 
happen at much higher collision energy, e.g. at CERN 
LHC? At what partonic scale the main processes will be 
stabilized? A plausible answer suggested in Ref. (2j]]3(J] is 
that high parton density will generate its own saturation 
scale, estimated for LHC to be about p 2 ~ 10 GeV 2 . 

If so, the instanton/sphaleron mechanism described 
above can no longer be important. The reason for that 
is extremely sharp dependence of the instanton effects 
on the scale involved, originating in semiclassical action 
exp(-S) ~ (A QC D/p) (11/3)Nc ~ {2/3)Nf - Therefore, if go- 
ing from RHIC to LHC we change p by factor 3, the 
sphaleron production is expected to drop by 3-4 orders 
of magnitude, becoming much less than its pQCD back- 
ground. 
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12 Although in QGP there are no quark condensates and one 
may think that all 6 't Hooft uuddss are produced, it is not 
nacessarily so since they could still be from the initial vacuum. 
Evaluation of probablities for each quark multiplicity we hope 
to report elsewhere. 
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